ICM 2002 • Vol. Ill • 1-3 



Quantum Resonances and 
Partial Differential Equations 



M. Zworski* 

Abstract 

Resonances, or scattering poles, are complex numbers which mathemati- 
cally describe meta-stable states: the real part of a resonance gives the rest 
energy, and its imaginary part, the rate of decay of a meta-stable state. This 
description emphasizes the quantum mechanical aspects of this concept but 
similar models appear in many branches of physics, chemistry and mathemat- 
ics, from molecular dynamics to automorphic forms. 

In this article we will will describe the recent progress in the study of 
resonances based on the theory of partial differential equations. 
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1. Introduction 

Eigenvalues of self-adjoint operators appear naturally in quantum mechanics 
and are in fact what we observe experimentally in many situations. To explain the 
need for the more subtle notion of quantum resonances we consider the following 
simple example. 

Let V(x) be a potential on a bounded interval, as shown in Fig^ a )- K £ 
denotes the classical momentum, then the classical energy is given by 

E=e+v{x), (i.i) 

and the motion of a classical particle can deduced from this equation by considering 
E as the Hamiltonian. The quantized Hamiltonian is given by 

P(h) = {hD x f + V{x) , £^hD Xl D x = -d x . (1.2) 
i 
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The operator P(h) considered as an unbounded operator on L 2 of a bounded in- 
terval (with, say, Dirichlet boundary conditions) has a discrete spectrum, which 
gets denser and denser as h — > 0, which corresponds to getting closer to classical 
mechanics: 

P(h)u{h) = E{h)u{h) , E(h)eR, J \u{h)\ 2 < oo , u(h) (±tt) = . (1.3) 

The eigenvalues, E(h), are what we (in principle) observe, and the square- integrable 
eigenfunctions, u{h), are the wave functions. 

Now, consider the same example but with a potential V(x) on R, as shown 
in Fig^b)- At the energy E, and inside the well created by the potential, the 
classical motion is the same as in the previous case. Hence, we expect that (at least 
for h very small) there should exist a quantum state corresponding to the classical 
one. It is clear in dimension one that for the potential shown in the figure, P (h) 
given by (|1.2(l has no square integrable eigenfunctions (|1.3J) . 




* n 



(a) (b) 

Figure 1: (a) A potential well on a finite interval, (b) A potential 
on the whole line the same classical picture but a very different 
quantum picture. 

A non-obvious remedy for this is to think of eigevalues (in the case of FigfTTkl 
as the poles of the resolvent of P(h): 

R(z,h) = (P(h)-z)- 1 . 

In the case of Fi stLb h the resolvent is not bounded on L 2 (R) for z > 0, which 
corresponds to free motion or scattering. However, under suitable assumptions on 
V(x) near infinity we have a meromorphic continuation of R(z, h) from Imz > to 
the lower half-plane: 

R(z,h) : Z4 mp (R) — iL(R-) • 

The complex poles of R(z, h) are the replacement for eigenvalues, are called reso- 
nances. For a recent presentation of these in the physical literature and for examples 
of current interest see for instance [7|. 

In this review, we will restrict ourselves to the case of 

P(h) = -h 2 A + V{x) , ieR", 
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where the potential V(x) satisfies some decay and analyticity assumptions near in- 
finity, or is simply compactly supported. Then our discussion above applies and 
resonances are defined as the poles of R(z, h) - see [§] for an attractive new presen- 
tation of the meromorphic continuation properties and references. In our convention 
the resonances are located in the lower half-plane. 

We should stress that most of the results hold for very general black box per- 
turbations of Sjostrand-Zworski t 22; which allow for a study of diverse problems 
without going into their specific natures - see [20] for definitions. Also, despite 
the fact that the motivation presented here came from molecular dynamics, sim- 
ilar issues arise in other settings, from automorpic scattering to electromagnetic 
(obstacle) scattering - see 0, and respectively. 

We suggest the the surveys ' GOI > > E3 > an d EH] for the review of earlier 
results, and we concentrate on the progress achieved in the last few years. For 
an account of work based on other methods and motivated by different physical 
phenomena we refer to [T3] . 

2. Trace formulae for resonances 

Trace formulae provide one of the most elegant descriptions of the classical- 
quantum correspondence. One side of a formula is given by a trace of a quantum 
object, typically derived from a quantum Hamiltonian, and the other side is de- 
scribed in terms of closed orbits of the corresponding classical Hamiltonian. A 
new general approach based on quantum monodromy operator which quantizes the 
Poincare map in a natural way was recently given in |24| . 

In general the spectral (or scattering) side of the formula is given in terms of the 
trace of a function of the quantum Hamiltonian f(P(h)). In the case of self-adjoint 
problems with discrete sets of eigenvalues the spectral theorem readily provides an 
expression for tr f(P(h)) but the problem becomes subtle for resonances. It has 
been studied by Lax-Phillips, Bardos-Guillot-Ralston, Melrose, Sjostrand-Zworski, 
Guillope-Zworski - see |34) and references given there. 

More recently Sjostrand |20 | .|21 j introduced local trace formulas for resonances, 
and that concept was developed further in |THj and 0. For P(h) given by (|1.2|l 
with V^a;) decaying sufficiently fast at infinity it can be stated as follows. 

Let n be an open, simply connected, pre-compact subset of {Re z > 0} C C, 
such that n R is connected. Suppose that / is holomorphic in a neighbourhood 
of O and that ip € (R) is equal to 1 in f2 n R and supported in a neighbourhood 
of O n R. Then, denoting the resonance set of P(h) by Res(P(h)), 

tr {(m(P(h)) - (fiP)(-h 2 A)) = ]T /(*) + E M (h) , 

zeRcs(P(h))nfi (2.1) 

\E M (h)\ < h- n C n ^max{\f(z)\ : z e Qi \ tt , Im z < 0} , 

where f2i is a neighbourhood of (see |H] and ^Hl for more precise versions). 
The basic upper bound on the number of resonances is given as follows: 



ft Res(P(/i)) n fi = 0(h- n ) , 



(2.2) 
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see ^21>EH>EI31>II21- It would consequently appear that the error term in (|2.1|) is of 

the same order as the sum over the resonances. However, by choosing the function 
/ so that it is small in fix \ O fl {Im z < 0}, the sum of f(z) can dominate the left 
hand side of Ij2.1|l . Doing that, Sjostrand has shown that an analytic singularity of 
E ^ \{x : V(x) > E}\ at Eo gives a lower bound in l|2.2|l for £1 a neighbourhood 
of Eq - see |23 an d references given there. 

Another application of local trace formulae techniques is in analysing reso- 
nances for bottles, that is perturbations of the Euclidean space in which the "size" 
of the perturbation may grow but which are connected to the euclidean infinity 
through a fixed "neck of the bottle" - see |21j.|18j. 



3. Breit-Wigner approximations 

In a scattering experiment the physical data is mathematically encoded in the 
scattering matrix, 5(A, h). It is in the behaviour of objects derived from the scat- 
tering matrix that we see physical manifestations of abstractly defined resonances. 
The classical and still central Breit-Wigner approximation provides this connec- 
tion. The most mathematically tractable case is provided by the scattering phase 
a(X,h) — logdet S(X, h)/(2iri), which is defined for T^'s with sufficient decay (oth- 
erwise relative phase shifts have to be used, see 0). When A is close to a resonance 
E — iT, the Breit-Wigner approximation says that 

^4 (A-i) 2 + r- (3-D 

that is, if r is small, cr(A) should change by approximately 1 as A crosses E. This 
has been justified rigorously in some situations in whch a given resonance, close to 
the real axis, is isolated. 

In view of (|2.2|l the number of resonances in fixed regions can be very large and 
for h small clouds of resonances need to be considered to obtain a correct form of 
the Breit-Wigner approximation. A formalism for that was introduced by Petkov- 
Zworski H5],|IE]) an d it was developed further by Bruneau-Petkov [3] and Bony- 
Sjostrand It is closely related to extending the trace formula 1|2.1|) to h dependent 
fi's. To formulate it, let us introduce ui(z,E) = (1/ir) J E (\lm z\/\z — X\ 2 )dX, Im z < 
0, the harmonic measure corresponding to the upper half-plane. Then, for non- 
critical A's (that is for A's for which £ 2 + V(x) = X => d (x ^(£ 2 + V{x)) ^ 0) 

cr(A + S, h) - a(X - S,h) = 0j(z,[X-5,X + 5]) + O(8)h- n , (3.2) 

zERe»(P(k)) 
\z-\\<h 

< 5 < h/C. The main point is that 5 can be made arbitrarily small and that we 
only need to include resonances close to A. If we do not assume that A is non-critical 
weaker, but still useful, results can be obtained from factorization of the scattering 
determinant |18| . or more generally, from the the analysis of the phase shift function 

m 

One of the consequencies of the development of the Breit-Wigner approxima- 
tions for clouds of resonances were new estimates on the number of resonances in 
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small regions, first in [T^], and then in greater generality in CD| , , CHI - If m place 
of (|2.2I) we had a Weyl law with a remainder then, as for eigenvalues, 

t( {z e Re$(P(h)) : \z - A| < 5} = 0{S)h-' 1 , Ch<6<l/C. (3.3) 

It turns out that despite the lack of the Weyl law we still have this estimate for 
non-critical A's. 

4. Resonance expansions of propagators 

In the case of discrete spectrum of a self-adjoint operator the propagator, 
exp(— itP(h)/h) can be expanded in terms of the eigenvalues. In fact, our under- 
standing of "state specificity" often comes from such "Fourier decompositions" into 
modes. For problems in which decay or escape to infinity are possible such expan- 
sions are still expected but just as in the case of trace formula? far from obvious. For 
non-trapping perturbations and in the context of the wave equation the expansions 
in terms of resonances were studied in late 60 's by Lax-Phillips and Vainberg (see 
the references in |30p. 

For trapping perturbations the expansions were investigated by Tang-Zworski 
|30| . and then by Burq-Zworski [5], Christiansen-Zworski 6 , Stefanov and 
most recently by Nakamura-Stefanov-Zworski |14| . 

Here, we will recall the general expansion given in let \ £ C-^°(Il n ), 
ip G C£°(0, oo), and let chsupp ip = [a, b). There exists < 6 < c(h) < 26 and L, so 
that we have 

xe -*«W^ (P(/l)) = £ X Res(e- lt '/ h R(.,h),z) X ij(P(h)) 

zeo(/i)nR CS (p) (4.1) 

+ L 2^ L 2(/i°°), for t > h- L , 

Q(h) = (a — c(h),b + c(h)) — i[0, 1/C) , and where Res(/(»), z) denotes the residue 
of a meromorphic family of operators, /, at z. 

The function c(h) depends on the distribution of resonances: roughly speak- 
ing we cannot "cut" through a dense cloud of resonances. Even in the very well 
understood case of the modular surface [HI Theorem 1] there is, currently at least, a 
need for some non-explicit grouping of terms. This is eliminated by the separation 
condition jSUl (4.4)] which however is hard to verify. 

The unpleasant feature of l|4.1|) is the need for very large times ~ h~ L and 
the presence of a non- universal parameter, c{h). The former is necessary in this 
formulation as one sees by considering the free case P(h) = —h 2 A. However an 
expansion valid for all times is possible in the case when a part of V constitutes a 
barrier separating the trapped set in {(x, £) : £ 2 +V(x) € supp ip} from infinity |14|. 
The example shown in FielTI b) is of that type. By the trapped set in S C T*R" 
we mean the set 

KnT, = {(x,£) e S : \exp(tH p )(z,g)\ -/-^ oo, t — > ±00} , (4.2) 



248 



M. Zworski 



where H p is the Hamilton vectorfield of p — £ 2 + V{x). One of the components 
comes from the work of Stefanov [23 on making estimates (|2.2|) more quantitative 
with constants related to the volume of the trapped set (for yet finer results of that 
type see Section. 6). 

5. Separation from the real axis 

One of the most striking applications of PDE techniques in the study of res- 
onances is the work of Burq (see 0] and references given there) on the separation 
of resonances from the real axis, estimates of the cut-off resolvent on the real axis, 
and the consequent estimates on the time decay of energy. 

For non-trapping perturbations we have the following estimate on the trun- 
cated meromorphically continued resolvent 

\\xR(z,h)x\\L^L* < C e ^°^ zl/h) , Imz>-MMogi, X eC™(R n ), 

(5.1) 

for any M, see for the absence of resonances (implicit in 1)5.1(1 ^. and ^1] for the 
(easily derived) estimate. 

Since the work of Stefanov- Vodev and Tang-Zworski (see [21] , EH] , |H3 an d 
references given there) we know that in many trapping situations 1 , there exist many 
resonances converging to the real axis 2 . The question then is how close could the 
resonances approach the real axis or, as it turns out equivalently, how big can the 
truncated resolvent be on the real axis. Heuristically, tunneling should prevent 
arbitrary closeness to the real axis, and the separation should be universally given 
by at least exp(—S/h), 

Tunneling of solutions is made quantitative in PDEs through Carleman es- 
timates used classically to show unique continuation of solutions of second order 
equations. Following the work of Robbiano, and Lebeau-Robbiano, Burq succeeded 
in applying Carleman estimates to a wide range of resonance problems - see |3] and 
references given there. In the setting described here his basic result says that if fl 
is as in (|2.1|) then 

3 5i,5 2 , z G Res(P(/i)) nfi Im z > -exp(-Si/h) , 
\\xR(z,h)x\\ L ^L*<exp(S 2 /h), zennR, x e C C °°(R"). 

In addition, improved estimates are possible if \ is assumed to have support outside 
the projection of the trapped set: 

\\ X R(z,h)x\\L^Li<C/h, zellnR, x e C™(TL n \ir(K)) . (5.3) 

That means that away from the interaction region tt(K) we have, on the real axis, 
the same estimates as in the non-trapping case (|5.1|l and that has immediate appli- 
cations in scattering theory . 

1 roughly speaking, whenever there exists an elliptic orbit of the //p-flow 

2 That in specific trapping situations there exist resonances whose distance to the real axis is 
of order exp(— S/h) is classical in physics with the most precise mathematical results given by 
Hclffcr-Sjostrand — see 0. Here we concentrate on very general existence results which guarantee 
lower bounds on the number of resonances. 
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Figure 2: Graph of a potential V(x), for which the classical flow 
of the Hamiltonian £ 2 + V(x) is hyperbolic for energies close to 
0.5. 



6. Resonances in chaotic scaterring 

Since the work of Sjostrand I19j on geometric upper bounds for the number of 
resonances, it has been expected that for chaotic scattering systems the density of 
resonances near the real axis can be approximately given by a power law with the 
power equal to half of the dimension of the trapped set (see below). Upper 

bounds in geometric situations have been obtained in |.'S2j and |35|. 

An example of a potential in V € Cj?°(R 2 ) for which the flow of H p , p = 
£ 2 + V(x), is hyperbolic (and hence scattering exhibits chaotic features) is shown 
in FiglrH 

A recent numerical study JUj for that potential indicates that the density of 
resonances satisfies a lower bound related to the dimension of the trapped set - see 
FigH^In a complicated semi-classical situation studied in |1U| . the dimension is a 
delicate concept and it may be that different notions of dimension have to be used 
for upper and lower bounds. That point is emphasized in |25j where numerical data 
for semi-classical zeta function for several convex obstacles is analyzed. 

In the case of convex co-compact hyperbolic quotients, X — T\H 2 , studied in 
|35| the situation is particularly simple as the quantum resonances coincide with 
the zeros of the zeta- function - see The notion of the dimension of the trapped 
set is also clear as it is given by 2(1 + 8). Here S = dimA(r) is the dimension of 
the limit set of T, that is the set of accumulation points of the elements of T (they 
are all hyperbolic), A(r) C OH 2 . 
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Figure 3: A plot of log(iV res ) as a function for — log(fi.) for the 
potential shown in in Fig. [JO The value of fi ranges from 0.025 
to 0.017. Triangles represent numerical data, circles least squares 
regression, and stars the slope predicted by the conjecture. 




Hence we expect that 

]T mr(s)~r 1+s , (6.1) 

|Im s\<r , Rc s>-C 

where mr(s) is the multiplicity of the zero of the zeta function of T at s. We also 
changed the traditional convention: here h 2 s(l — s) = z in the notation of previous 
sections, and h(Re s — 1/2), Mm s correspond to Im z, Re z respectively. 

An upper bound of this form was established in but a simpler method 
giving improved upper bounds has been recently presented in |37| . The new method 
is based on zeta function techniques. What we obtain is a bound in the case of 
convex co-compact Schottky groups 

^{m r (s) : r<Ims<r + l, Re s > -C } < Cir s , (6.2) 

where 5 = dim A(r). This improved estimate is a "fractal" version of l|3.3|) : 1 + 5 
now plays the role of n. It is clear that the method works in greater generality and 
implementing it is part of an ongoing project. 
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